Signatures of Prelocalized States in Classically Chaotic Systems 
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to TRS while a direct optimal fluctuation (DOF) 

method predicts a symmetry independent result M . Re- 
cent numerical calculations j9) on the Anderson model 
seem to support the latter theory. This controversy, was 
an additional motivation for the present work. 

In the present article, we numerically study the distri- 
bution of intensities of the Floquet- states of the kicked 
rotor (KR) on a torus |l6| and its 2d generalization jL7) . 
Our system is defined by the time-dependent Hamilto- 



H = H a + kVY,5(t- mT) , Ho({&}) = £ + 7 i) 
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We investigate the statistics of eigenfunction intensities 'PQV'I 2 ) in dynamical systems with clas- 
sical chaotic diffusion. Our results contradict some recent theoretical considerations which challenge 
the applicability of field theoretical predictions, derived in a different framework for diffusive dis- 
ordered samples. For two-dimensional systems, the tails of 7'(|V'| 2 ) contradict the results of the 
optimal fluctuation method, but agree very well with the predictions of the non-linear a— model. 
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The statistical properties of wavefunction intensities 
have sparked a great deal of research activity in recent 
years. These studies are not only relevant for mesoscopic 
physics ]l|-[TrH , but also for understanding phenomena in 
areas of physics, ranging from nuclear ]Tl[ | and atomic 
p2[ to microwave physics and optics Hj. Experi- 
mentally, using microwave cavity technics it is possible to 
probe the microscopic structure of electromagnetic wave 
amplitudes in chaotic or disordered cavities |l3|]. Re- 
cently, the interest in this problem was renewed when 
new effective held theoretic techniques were developed 
for the study of the distribution of eigenfunction inten- 
sities 'P(|'0| 2 ) of random Hamiltonians. As the disorder 
increases, these results predict that, the eigenfunctions 
become increasingly non-uniform, leading to an enhanced 
probability of finding anomalously large eigenfunction in- 
tensities in comparison with the random matrix theory 
(RMT) prediction. Thus, the notion of prelocalized states 
has been introduced Jt)-^ to explain the appearance of 
long tails in the distributions of the conductance and 
other physical observables [jD. 

Up to now all theoretical predictions and numer- 
ical calculations ||[l0) apply to disordered systems and 
are based on an ensemble averaging over disorder real- 
izations. Their validity, however, for a quantum dynam- 
ical system (with a well defined classical limit) that be- 
haves diffusively is not evident. Furthermore, based on 
an argument put forward in M (see also JL5||), the far 
tail of V(\ip\ 2 ) is due to rare realizations of the disorder 
potential, and therefore requires an exponentially large 
number of eigenfunctions, which can only be accounted 
by disorder averaging. Here instead we study the statis- 
tical properties of eigenfunctions in a dynamical model 
without introducing any ensemble averaging. Our main 
conclusion is that in a generic dynamical system with 
classical diffusion, "Pd^l 2 ) is described quite well by the 
nonlinear cr— model (NLSM). We point out here that be- 
tween the various theoretical works there is a consider- 
able disagreement about the parameters that control the 
shape of 'P(|V'| 2 ) an d their dependence on time-reversal 
symmetry (TRS). More specifically, the NLSM suggests 
that the tail of P(\ip\ 2 ) in two dimensions (2c?) is sensitive 



V({9i}) = cos(0i) cos(6> 2 ) cos(a) + - sin(26>i) cos^a) sin(a) 

where £, denotes the angular momentum and 0i the con- 
jugate angle of one rotor. The kick period is T, k is the 
kicking strength, while Tj is a constant inversely propor- 
tional to the moment of inertia of the rotor. The standard 
KR corresponds to d — 1 (with 82 = 0) whereas for d = 2 
we have a two-dimensional generalization. The parame- 
ter a breaks TRS |r|, the parameters ji are irrational 
numbers whose meaning will be explained below. The 
Hamiltonian ([!]) describes a system which is kicked pe- 
riodically in time and is integrable in the absence of the 
kicking potential. The motion generated by (|l|) is classi- 
cally chaotic and for a sufficiently strong kicking strength 
k there is diffusion in momentum space with diffusion co- 
efficient D = limt^oc < £ 2 (t) > /t ~ fc 2 /2T (within the 
random phase approximation) 16|l^ |. 

If the Ci are taken mod(27rmi/Tri) where rtii are in- 
tegers, Eq. (|l|) defines a dynamical system on a torus. 
The quantum mechanics of this system is described by a 
finite-dimensional time evolution operator for one period 



U = exp(-iH ({d})T) exp{-iV({ei})) 



(2) 



where we put Ti — 1 . Upon quantization, additional sym- 
metries associated with the discreteness of the momen- 
tum show up, which can be destroyed by introducing ir- 
rational values for the parameters 7j's. The most striking 
consequence of quantization is the suppression of classical 
diffusion in momentum space due to quantum dynamical 
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localization [ Jl€3|Jl 7| ] . We introduce the eigenstate compo- 
nents Vl/fc(n) of the Floquet operator in the momentum 
representation by 



(3) f 



The quantities Wfc are known as quasi-energies, and their 
density is p = T/2n. The corresponding mean quasi- 
energy spacing is A = l/(pL d ), where L is the linear size 
of the system. The Heisenberg time is tn = 2ir/ A while 
to = L 2 /D is the diffusion time (Thouless time). Now 
one can formally define a dimensionless conductance as 
g = tn/to = DkL d ~ 2 where Dk = TD is the diffusion 
coefficient measured in the number of kicks. Four length 
scales are important here: the wavelength A, the mean 
free path Im, the linear extent of the system L, and the 
localization length £. According to Refs. §-0] the field 
theoretical predictions are derived under the conditions 



A < l M < L < £. 



(4) 



The first condition ensures that transport between scat- 
tering events may be treated semiclassically. This limit 
can be achieved for our system ([!]) when k — > oo, T — > 
while the classical parameter K — kT remains constant. 
When Im <C L as long as the motion is not localized 
(i.e. L <C Ct it is diffusive, since a particle scatters many 
times before it can traverse the system. The resulting 
mean free path for our system is Im ~ y/~D~k while 
the localization length for d — 1 is £ ~ -Dfe/2 and for 
d=2is£~ lMe flt/2 @@. 

Here we calculate the distribution function V(t = 
L d \^fk(n)\ 2 ) jl9| by using a direct diagonalization of the 
Floquet operator (||). The TRS is broken entirely for 
a = 5.749. In order to test the issue of dynamical cor- 
relations, we randomize the phases of the kinetic term 
of the evolution operator (||) and calculate the result- 
ing V(t). This model will be referred to as Random 
Phase KR (RPKR). Since all our eigenfunctions have the 
same statistical properties (in contrast to the Anderson 
cases where one should pick up only eigenfunctions hav- 
ing eigenenergies within a small energy interval [p|Jlo|) 
we make use of all of them in our statistical analysis. 
The classical parameter K is large enough in all cases 
to exclude the existence of any stability islands in phase 
space. The classical diffusion coefficient Dk is calculated 
numerically by iterating the classical map obtained from 
(Q) . Below we present our numerical results and compare 
them to the predictions of Refs. 

Id Kicked Rotor. It was shown in jla], that the ef- 
fective field theory describing the semiclassical physics 
of the system is precisely the NLSM for quasi-one di- 
mensional (Id) metallic wires. Such a mapping however, 
requires an averaging over an ensemble of rotors having 
the same classical limit. We point out again that in the 
calculations below we do not adopt such an averaging 
procedure. 
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FIG. 1. Corrections to the distribution intensities &Vp{t) for 
the kicked rotator model i.e. Eq. (Q) for d = 1. The system size is 
L = 1024, (o) p = 1, (O) ft = 2. The solid (dashed) lines are the 
best fit of (|) for fi = 1(2) to the numerical data: (a) D k w 1800 
and (b) D k fa 3150 ; (c) Shows the extracted diffusion propagator 
K@ vs. L/D k . 

The NLSM for quasi-ld systems can be solved ex- 
actly for the distribution function "Pp{t), using a transfer 
matrix approach l^,^. I n t ne ballistic regime (where 
g -> oo) RMT is applicable and one finds V%£fdt) = 

exp(-t/2)/V2ni and 7>j^(f) = exp(-i) @. Here (3 
denotes the corresponding Dyson ensemble ((3 = 1(2) for 
preserved (broken) TRS). As localization increases, the 
deviations from the RMT results of the body and the 
tails of the distribution Vp(t) become noticeable and can 
be parameterized by a single parameter which is the di- 
mensionless conductance g = Dk/ L. 

For t < yjDk/L, according to all studies |3|-|3| V(t) 
is just the RMT result with polynomial corrections in 
powers of L/D k , i.e. Vp{t) = Vf MT {t)[l + 5Vp(t)}. The 
leading term of this expansion is given by 



SPp(t) 



( 3/4-3i/2 + f 2 /4 for p= 1 
K \l-2t + t 2 /2 for/? = 2 



(5) 



where k ~ 1/g is the fd diffusion propagator, which is 
identical for (3=1 and (3 = 2 since it is a classical quan- 
tity. 

In Fig. la,b we report our numerical results for SVp(t) 
for two representative values of Dk- One can clearly see 
that the agreement with the theoretical prediction (^) 
becomes better as Dk increases. This is due to the fact 
that by increasing Dk we are approaching the semiclas- 
sical region and therefore Eqs. (0) are better satisfied. 
At the same time higher order corrections in SVpif) be- 
come negligible with respect to the leading term given 
by Eq. (j^). The resulting k\ and K2 obtained by the 
best fit of our data to Eq. rt|) are found to be equal 
and in excellent agreement with the theoretical value (see 
Fig. Ic). We therefore conclude, that in a generic dynam- 
ical system, the only parameter that controls the shape 



2 



of the deviations 5Vp(t) is the classical diffusion propa- 
gator. Moreover, our results are in excellent agreement 
with the recent NLSM predictions derived in the frame- 
work of diffusive disordered systems. Finally in Fig. lc 
we also report the outcome of the RPKR model. The 
results remain essentially the same indicating that Vp (t) 
for quasi-ld systems are insensitive to dynamical corre- 
lations. 
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FIG. 2. (a) Tails of the distribution P@(t > Dk/L) for the 
model (0) for d = 1 with L = 1024, D fc ~ 2625 and for /3 = 1 
(o) and = 2 (O). The solid (dashed) lines are the best fit of (§) 
for /3 = 1(2) to our data; (b) Coefficients C 5 vs. D k /L. The 
solid (dashed) lines are the best fits to Cp = ApW Dk/L + Bp for 
= 1(2). 

The tail of the distribution [t > Dk/L) deviates 
strongly from the RMT prediction and has a stretched 
exponential form |j] || 



Vp(t) ~ ^ exp(-2CpV~t), - P^/D k /L (6) 

where ^ is a symmetry dependent constant. Our nu- 
merical results agree nicely with Eq. (^). In Fig. 2a 
we present an example of Vp(t). By fitting our data 
to Eq. the coefficients C\^Ci can be extracted. In 
Fig. 2b we report the extracted stretched exponential co- 
efficients Cp from the best fit of (JsJ) as a function of the 
square root of the dimensionless conductance g — Dk/L. 
A nice linear behavior is observed. The best linear fit 
Cp = Apy/D k /L + Bp yields, Ap =1 = 0.41 ± 0.05 and 
Ap =2 = 0.82 ± 0.05. The resulting ratio R = A 2 /Ai = 2 
is in excellent agreement with the theoretical prediction 
(H). We have also calculated the stretched exponential 
coefficients Cp for the RPKR model. The results for var- 
ious Dk values are summarized in Fig. 2b and show a 
nice agreement with the results obtained from the real 
Hamiltonian. 

2d Kicked Rotor. According to Ref . || , corrections to 
the body of J>R- MT are still given by Eq. (||) , but now k 
is the 2d diffusion propagator. 

Figures 3a, b show corrections to , p^- MT for g = Dk ^> 
1 for two representative values of Dk- We find again 



that the form of the deviations are very well described 
by Eq. (||) and the agreement becomes better for larger 
values of the diffusion constant. In Fig. 3c we summa- 
rize our results for various Dk values. The extracted Kp 
values are obtained by the best fit of the data to Eq. (||). 
Again we find that Kp depends linearly on l/Dk- How- 
ever, contrary to the ld-KR, here K\ and K 2 , are differ- 
ent. Moreover the best fit with Kp = ApD^ 1 + Bp yields 
Ap =1 = 5.44 ± 0.03 and Ap =2 = 10.84 ± 0.04 indicating 
that the ratio R = A 2 /A\ is close to 2, a value that could 
be explained on the basis of ballistic effects . Tak- 
ing the latter into account leads to an additional term in 
the classical propagator Kp = Kdiff + 2 K baii- The first 
term is the one discussed previously and is associated 
with long trajectories which are of diffusive nature while 
the latter one is associated with short ballistic trajecto- 
ries which are self-tracing ||. Thus, when Kdiff ^ K baii 
we get R = 2. The calculation with the RPKR model 
shows, however, that the corresponding ratio is R ~ 1 
in agreement with the theoretical prediction for disor- 
dered systems with a pure diffusion. This indicates that 
dynamical correlations can be important in the 2d case. 
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FIG. 3. Corrections to the distribution intensities SVp(t) for 
the kicked rotator model ([l]) for d = 2. The system size is L = 90, 
(o) = 1, (O) /3 = 2. The solid (dashed) lines arc the best fit 
of (^|) for /3 = 1(2) to the numerical data: (a) Dk ~ 34 and (b) 
Dk ta) 53 ; (c) Fit parameters Kp vs. D^ 1 . The solid (dashed) lines 
are the best fits to Kp = ApD~ x + Bp for f} = 1(2). 

For the tails of the distributions, the result of the 
NLSM within a saddle-point approximation is 



Vp{t) ~ exp[-<7£(lnt) 2 ], C% 



(3tt 2 P D 
2 ln(L/0' 



(7) 



Note that the decay in the tails of Eq. (m) depends on 
0, as in the l<i-KR case (see Eq. ([}])). Recently, a DOF 
method was used to calculate the tails of Vp(t) JpJ. It 
was found that the tails are still given by Eq. (ffl) but 
with a log-normal coefficient C which is independent of 
the parameter j3 : 
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DOF 



7T 2 p 



D 



ln(L/A)' 



(8) 
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FIG. 4. (a) Tails of the distribution Vp =1 {t > D k ) for the 
model (Q) for d = 2 and Dj. ~ 35. The system size is L = 80, (o) 
f3 = 1, (O) (3 = 2. The solid (dashed) lines arc the best fit of (|) for 
/3 = 1(2) to the numerical data; (b) Fitted log- normal coefficients 
Cp versus the classical diffusion coefficient D k . The solid (dashed) 
lines are the best fits to Cp = ApD k + Bp for /3 = 1(2). 

Figure 4a shows a representative case of Vp = i{t > 
Dk). The tails show a log-normal behavior predicted by 
Eq. ([?]). In Fig. 4b we report the log-normal coefficients 
C/3 extracted from the best fit to our numerical data, 
versus the classical diffusion coefficient. A pronounced 
linear behavior is observed in agreement with both theo- 
ries. However one clearly sees that C\ differs from C2 in 
contrast to the DOF prediction (|^) and to recent numer- 
ical calculations done for the 2d Anderson model j^] . We 
point out here that in |9] the authors were not able to go 
to large enough values of conductance g (in comparison 
to our study) where the theory can really be tested. In 
contrast, the NLSM predicts a value of 2 for the ratio 
R = Cf/Cf . We note that Cp is only the leading term 
in Dk- In order to calculate this ratio, we performed a 
fit to our data with Cp = ApDk + Bp. The resulting ra- 
tio was found to be R — A2/A1 = 1.97 ± 0.03 in perfect 
agreement with the NLSM predictions. Finally in Fig. 4b 
we also present our results for the RPKR model (using 
the same data as the one in Fig. 3d). Again we found 
that the ratio R = 1.96 ± 0.03 w 2. Thus V{t > D k ) 
depends on TRS and is described by the NLSM. 

In summary, we have performed a detailed numerical 
analysis of the statistical properties of the wavefunction 
intensities V(t) of the standard KR on a torus and its 
2d generalization. Based on these results, we concluded 
that the distribution V{t) of generic quantum dynam- 
ical systems with diffusive classical limit is affected by 
the existence of prelocalized states. The deviations from 
RMT are well described by field theoretical methods de- 
veloped for disordered systems. In particular, in a clari- 
fying way we have resolved the controversy between DOF 



and NLSM by demonstrating that the dependence of the 
tails of Vpit) on TRS is described correctly by the latter 
theoretical approach. 

We acknowledge useful discussions with L. Kaplan. 
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